Higher-order supersymmetric quantum mechanics 



David J. Fernandez C. and Nicolas Fernandez-Garcfa 
Departamento de Ffsica, CINVESTAV 
AP 14-740, 07000 Mexico DF, Mexico 

in 
o 
o 

CsJ 1 Abstract 

We review the higher-order supersymmetric quantum mechanics (H-SUSY QM), 
\ which involves differential intertwining operators of order greater than one. The it- 

■ erations of first-order SUSY transformations are used to derive in a simple way the 

higher-order case. The second order technique is addressed directly, and through this 
approach unexpected possibilities for designing spectra are uncovered. The formalism 
is applied to the harmonic oscillator: the corresponding H-SUSY partner Hamiltonians 
are ruled by polynomial Heisenberg algebras which allow a straight construction of the 
^ ! coherent states. 
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1 Introduction 



The number of exactly solvable problems in non-relativistic quantum mechanics is small, 
and most of them can be dealt with the factorization method. This technique, introduced 
long ago by Schrodinger [1,2], was analyzed in depth by Infeld and Hull [3], who made an 
exhaustive classification of factorizable potentials. Later on, Witten noticed the possibility 
of arranging the Schrodinger's Hamiltonians into isospectral pairs (supersymmetric part- 
ners) [4]. The resulting supersymmetric quantum mechanics (SUSY QM) revived the study 
of exactly solvable Hamiltonians (see e.g. [5]). An additional step was Mielnik's factoriza- 
tion through which the general SUSY partner for the oscillator for a certain factorization 
energy was found [6]; this technique was immediately applied to the hydrogen potential [7]. 
Meanwhile Nieto [8], Andrianov et al [9] and Sukumar [10-14] put the method on its natural 
background discovering the links between SUSY, factorization and Darboux algorithm, caus- 
ing then a renaissance of the related algebraic methods [15-55] (for a recent review see [56] 
and references therein). 

These procedures can be recovered from a general setting in which a first-order dif- 
ferential operator intertwines two Hamiltonians [57,58]. This so-called first-order inter- 
twining technique suggests further generalizations: the most obvious one involves a /c-th 
order differential intertwining operator and gives place to the higher-order supersymmetric 
quantum mechanics (H-SUSY QM) [59-77]. It can be achieved by iterations of first order 
SUSY transformations [64]. An alternative to deal with the problem consists in looking 
for the k-th order operator directly, expressing the intertwiner as a sum of the k + 1 terms 
gi(x)d l /dx l ,i = 0, . . . , k, and solving the system of equations resulting from the intertwin- 
ing relationship for the gi(x)'s. By assuming that the initial Hamiltonian is solvable, new 
solvable ones and their eigenstates can be generated through the two previous formulations. 

Concerning applications, the procedure works successfully to generate /c-parametric fam- 
ilies of potentials almost isospectral to the harmonic oscillator [64,67], the radial hydrogen- 
like potentials [65,66], and in the free particle case [69] (for a collection of recent papers on 
SUSY QM see [78]). Specially interesting are the SUSY partners of the oscillator because 
through them some connections with other important subjects of mathematical physics can 
be established. 

In the first place let us notice the existence of differential annihilation and creation 
operators of order greater than one for the SUSY partners of the oscillator. These operators 
provide natural realizations of the so called polynomial Heisenberg algebras [67,75,79-83]. It 
is interesting as well that these non-linear algebras admit partial linearizations [84-86], i.e., 
through appropriate modifications on the ladder operators the standard Heisenberg algebra 
is recovered on the subspace spanned by the eigenf unctions intertwined with the physical 
eigenstates of the oscillator. 

In the second place let us mention the construction of coherent states (CS) for potentials 
generated through the intertwining technique [67,82,84-94]. In this direction Fukui and 
Aizawa have derived some CS for the 'shape invariant' potentials [87], a particular class of 
solvable potentials generated long ago by Infeld and Hull through the factorization method 
[3]. The first work involving CS for the simplest non-trivial family of potentials isospectral 
to the oscillator (Abraham-Moses [6,95]) was done in 1994 [84] (see also [85,86]). Later on 
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various developments have appeared, e.g., some authors have constructed the CS for a class 
of anharmonic oscillators which spectra consist of a part isospectral to the oscillator plus one 
lower energy at a multiple of the spacing between the oscillator levels [91]. Furthermore, CS 
have been derived [82] for families of SUSY partner potentials having an arbitrary eigenvalue 
below the ground state energy of the oscillator [11,43,64]. Recently, it has been implemented 
as well the CS construction for the general H-SUSY partners of the oscillator [67]. 

In these lecture notes we are going to address the subjects mentioned above. We will 
start with the iterative H-SUSY QM to generate families of potentials (almost) isospectral 
to the initial one. Then, we will formulate the same problem directly, by assuming that the 
intertwining operator is of k-th order. Due to the difficulty involved in this problem, we just 
will illustrate the technique through the second-order supersymmetric quantum mechanics. 
Then we will apply these generation procedures to the harmonic oscillator. It will be shown 
that the intrinsic symmetry of the H-SUSY partners of the oscillator is generated by a pair 
of differential annihilation and creation operators of order greater than one, giving place 
to explicit representations of the polynomial Heisenberg algebras. We will show as well an 
interesting linearization procedure for the nonlinear algebras characteristic of the H-SUSY 
partners of the oscillator. We will perform then the coherent states construction in the non- 
linear as well as in the linear cases. We will conclude these notes with some discussion on 
the future of SUSY QM. 
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2 Higher-order supersymmetric quantum mechanics: 
iterative approach 

Let us start with two Schrodinger Hamiltonians 

and suppose the existence of a first order differential operator A\ intertwining them 

H X A\ = A\H (2) 

where 

A = As(-4- + a i( x )) ( 3 ) 



V2 \ dx 

the superpotential ai(x) being a real function to be determined. In order to 'pass' the 
differential operator d/dx and its powers to the right in the terms arising from (j2J, we use 
the following operator relationships 

ax ax 

d 2 d 2 d 

aar aar ax 

Then, it is straightforward to show that 



Id 3 aci d 2 . d 



Due to (j2J) we make equal the coefficient of the same powers of d/dx of the two previous 
equations to obtain 

Vx = V Q - a[ (4) 

aiVx - ^ = aiF - Vq' (5) 

Substituting the expression for V\ of (jU in © and integrating the result we get: 

«; + a? = 2{V - e) (6) 

From now on we are going to express explicitly the dependence of the superpotential in 
terms of the factorization energy e in the way ai(x,e). If a function u^°'(x) such that 
ai(x,e) = /u^ is used, we have that (|4J6|) become: 

A u (o)" + VqU (o) = m (o) (7) 



i.e., is a solution (not necessarily physical) of the initial stationary Schrodinger equation 
associated to e. 

Let us notice that ()4I6|) guarantee that H and Hi become factorized: 

H = A X A\ + e H 1 = A\A X + e (9) 

where 

Al = T2(^ + ai(x ' f) ) (10) 

is the operator adjoint to A\. 

Suppose that Vq{x) is a solvable potential with eigenfunctions ipn\x) and eigenvalues 
E n , n = 0, 1, . . . Furthermore, let us assume that we have found a solution ai(x, 6i) (i4° (x)) 
to the Riccati equation © (Schrodinger equation (J2J)) for a given value of the factorization 
energy e = ei < E'o, where E is the ground state energy of H . Thus, the V\{x) given in (jlj) 
is a completely determined solvable potential with normalized eigenfunctions 

^)( a .) {xe -J7«i(».«)* = -^_ (ii) 

Itj (x) 

^ 1} (x) = ^M (12) 

and eigenvalues {ei, E n , n — 0, 1, . . . }. Let us remark that the restriction t\ < E is crucial 
to avoid the existence of singularities in ai(x, ei), in Vi and also in the ^ of (fTTfT^ . 
Indeed, if €\ would be greater than E , the transformation function would have always 
zeros in the initial x-domain and thus «i(x, ei) would have singularities at those points. If 
however e\ < Eq, then u^\x) will have at most one zero; by exploring the two-dimensional 
subspace of solutions associated to e\ it is possible to find a subset of nodeless solutions (see 
e.g. [10,11]). By simplicity, we shall assume that the factorization energy used to generate 
any new Hamiltonian through the first-order SUSY is below the ground state energy of the 
initial Hamiltonian. We shall suppose as well that, for fixed e±, the arbitrary parameter of a 
general Riccati solution © has been adjusted in order to avoid the singularities in the a's. 

Now we iterate the previous technique, taking V\(x) as the given solvable potential to 
generate a new one V^(x) through an intertwining operator A\ and a different factorization 
energy €2, with e 2 < ex. The corresponding intertwining relationship, H 2 A\ = A 2 Hi, leads 
to equations similar to (|4I6|) : 

a' 2 (x, e 2 ) + al(x, e 2 ) = 2[Vi(x) - e 2 ] (13) 
V 2 (x) = V 1 (x)-a' 2 (x,e 2 ) (14) 

In terms of transformation functions u 2 \x) such that a 2 (x, e 2 ) = ju^ 1 we have: 

A u m" + VlU w = e2U w ( i5) 



V, = W-[\\ (16) 
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An important result to be proved now is that the solutions to (fT3*j) can be algebraically 
determined through solutions of the initial Riccati equation (JSJ for the factorization energies 
ei and e 2 [64-69]. We stick to [96] because there the calculations are transparent. To find 
the corresponding formula [69], we know that the two initial Riccati solutions satisfy: 

a[(x, ex) + a\(x, ex) = 2[V (x) - ex] 

a[(x, e 2 ) + a\{x, e 2 ) = 2[V (x) - e 2 ] (17) 



At the Schrodinger level we have that 

HqUx(x) = exuf\x) 
H Q u 2 \x) = e 2 u 2 \x) 

where 

uf\x) oc e ^ ai ^ ei]dy 
u [ 2 \x) oc e ^ ai{y ^ )dy 



Let us remember that u\ (x) is used to implement the first SUSY transformation and that 
the eigenfunction of Hx associated to ex is given by (jllj) . On the other hand, the eigenfunction 
of Hx associated to e 2 is given by: 



ytrl.M „ (°)\ 

(1) a\ (0) (0)' / n (0) ,u 2 ) 

u 2 oc A[u 2 oc — u 2 + ax{x, ex)u 2 oc ^ (19) 



and since that 



(0)' / \ (o) 

u 2 = cxi{x, e 2 )w 2 



we get: 

u 2 oc [«i(x, ex) — ctx(x, e 2 )] u 2 (20) 

In order to implement the second SUSY transformation, we express in the standard way 
in terms of the corresponding superpotential: 

u { 2 \x) oc e ^ a ^ €2)dy (21) 

By plugging (j2~Tj) in (J2U|l we arrive at 

e f *a 2 (y,e 2 )dy ^ ^ _ a ^ ^ ^(0) ^ 

By taking the logarithm of both sides of the previous equation: 

/ a 2 {y, e 2 )dy = In u 2 + In [«i(x, ei) — ax(x, e 2 )] + constant (23) 
Jo 



Deriving this expression with respect to x: 



ot 2 (x, e 2 ) = e 2 ) + 



a[{x,€i) - a[(x, e 2 ) 



a x [x,ei) - a 1 {x,e 2 ; 
By using the initial Riccati equations (|17jl we obtain: 

a'^x^i) - a[(x,e 2 ) , N / , 2(ei - e 2 ) 



(24) 



oti{x,ei) - a 1 (x,e 2 ) 



—a>i{x, e 2 ) - a>i{x, t\) - 



a 1 (x,e 1 - a 1 lx,e 2 ) 



Therefore 



a2(x,e 2 = -ai(x,eij - 



2(e a -e 2 ) 



(25) 



(26) 



oci\x,ei) - ai{x,e 2 ) 

This formula expresses the solution to (fT3Jl with Vi(x) = V (x) — a[(x, e x ) in form of a finite 
difference formula involving two solutions ai(x, ei), «i(x, e 2 ) of the Riccati equation (jOJ) for 
the factorization energies ei, e 2 (see also [64]). Notice that a similar formula has been used 
by Adler in order to discuss the Backlund transformations of the Painleve equations [80]. 
The potential V 2 (x) reads: 



V 2 (x) = V x {x) - a' 2 (x, e 2 ) = V (x) + 

The eigenfunctions associated to H 2 are given by: 

ip£\x) oc e -Jo a *(y^)dy = . 



2(ei-e 2 ) 



ai(x,ei) - ai(x, e 2 ) 



(27) 



u 



(i) 



(2f 



V^n - e 2 v/(£„-ei)(£„-e 2 ) 



The corresponding eigenvalues are {e 2 , ei, _E n , n = 0, 1, . . . }. 

We can continue the iteration process as many times as solutions for different values q 
to the initial Riccati equation © we have. Suppose that we know k of these, {ax(x, i — 
1, . . . , k, Cj+i < 6i}, and we iterate the process k times. Hence, a new solvable Hamiltonian 
Hk will be gotten whose potential reads: 



V k (x) = V k -i(x) - a' k (x, e k ) = V (x) - a-(x, e;) 



(29) 



i=i 



where a«(x, e«) is given by a recursive finite difference formula generalizing (j2EJ): 

2(e* - e i+ i) . 



k-1 



(30) 
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The eigenfunctions are given by: 



1p£\x) OC e ~ Jo a k(y,e k )dy 

0) = 



(31) 



V (ei-e2)...(ei-£fc) 
v /(£„-e 1 )...(£: n -e fe ) 



The corresponding eigenvalues are {q, i — k, . . . , 1, n = 0, 1, . . . }. 

In order to have the scheme complete, let us remember how are intertwined the H^s: 

H i A\ = A\H^ 1 i = l,...,k (32) 

Thus, departing from H we have generated a chain of factorized Hamiltonians: 

H t = A\Ai + e { = A i+1 A\ +1 + e l+1 i = 1, . . . , A; - 1 (33) 
H k = A\A k + e k (34) 

where the end potential V k (x) can be recursively determined by means of (|29|3()j) . Thus, we 
require the mentioned k solutions ai(x, e*), i — 1, . . . , k to the initial Riccati equation (0), 
which means to have k non-equivalent factorizations of the Hamiltonian H : 

H ^^(^ + a 1 (x,ei)j(-^ + a 1 (x,ei)j+e i i = l,...,k (35) 

Let us notice that there is a fc-th order differential operator, B\ = A\ . . . A±, intertwining 
the initial Hq and final Hamiltonians H k : 

E k B\ = BlH (36) 

From equations (|ST|) we get: 

Bt^ = ^ {En _ ei) . . . {En _ €k) (37) 

while from the adjoint to (jHEj) it turns out that: 



B k ^ n k) = ^{E n -e 1 )...{E n -e k )^ (38) 

These equations lead to the higher-order supersymmetric quantum mechanics H-SUSY QM 
[59-69]. In this treatment the standard SUSY algebra with two generators [4] 

[Qi, H ss ] = {Q t , Qj} = 5 i:j H SB i,j = 1, 2 (39) 
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is realized with the aid of B k and B\: 



(40) 



(41) 



where Q x = (Qt + Q)/^2, Q 2 = i(Q - Q j )/V2. Due to 

4s fc = A\...A\A l ...A k = A\...A\{H 1 -e 1 )A 2 ...A k 

= A[...Al(H 2 - ei)(ff 2 - e 2 )A 3 . . . A fe = (if fc - ei) . . . (H k - e k ) 

B k B\ = A l ...A k A\...A\ = A l ...A k _ l {H k _ l -e k )A\_ 1 ...A\ 

= (H -e 1 )...(H -e k ) 



it turns out that the SUSY generator H ss is a fc-th order polynomial 



of the physical Hamiltonian H% involving the intertwined Hamiltonians H and H k 



In particular, if k — 1 we will get the standard SUSY QM for which H ss = (H% — e x ), i.e., 
H ss is linear in Hf. If k = 2 we will get the quadratic superalgebra, or SUSUSY QM due to 



H ss = - Cl ) . . . (HI - e k ) 



(42) 




(43) 



H ss = (Hi - ei)(if? - e 2 ) [59,62,63]. 
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3 Higher-order supersymmetric quantum mechanics: 
direct approach 

In the direct procedure to the H-SUSY QM it is supposed from the very beginning that the 
intertwining operator in is of k-th order: 



Bl = 2~ k ' 2 



d k d k ~ l d 

l) fc 7T + 9k-iix)-—— + ■ ■ ■ + gi(x)— + g {x) 



' dx k i " t ±v ' dx k 1 ' aJA ' dx 



(44) 



where the real functions {gi(x),i = 0, . . . , k — 1} can in principle be determined through an 
approach similar to the one followed in the first-order case. Let us notice that the formulae 
(j36H43|) are still valid in this case, but the procedure to derive them is long and involved. 
Here we present just the simplest case with k = 2. This will illustrate clearly the advantages 
to design spectra supplied by the direct procedure compared with the iterative one. 

3.1 Second-order supersymmetric quantum mechanics 

The second-order supersymmetric quantum mechanics [59-62,77,96] arises when considering 
a second order intertwining operator such that 



H 2 B\ 


= B\H { 
1 d 2 


Hi = 




~2~ch? 


B\ = 


1 / d 2 


2 [dx 2 ' 



(45) 



+ Vi(x) z = 0,2 



dx 



The calculation of the left hand side of (|4*3j) leads to: 



t Id 4 rj d 3 f , 7 „\ d 2 



2H ^ = -2^ + i^H"'-i + ^^ 



+ [^_y_,V,)| + 7 14-r (46) 



The corresponding right hand side provides: 



+ (2\/J - rjV ) ^ + V£ - + 7 Vo (47) 



+ (y 



52 "° ~2dx~^2dx~^V°~2Jdx^ 

d 

dx 



The coefficients of the same powers of d/dx in ()46I47|) must be equal, and thus we arrive at: 

V 2 = V - i (48) 
!f - y - V V 2 = 2V ' - 77V0 (49) 
lV 2 - dL = V £ - V V> + 7V0 (50) 
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By substituting in and solving for 7' we get: 



i = \ + W - 2^0' (51) 

Integrating this equation with respect to x: 

1 2 

1 = 1 L + r L-2V, + d (52) 

where d is a real constant. By plugging in (J5Uj) and solving for 7" we have: 

7" = 2 (rjV^ - V£ - jrf) (53) 
Taking the derivative of (JSTj) with respect to x we get: 

7" = ^ + W + V 2 - 2y " (54) 



If we substitute ()52|) and (J53)l in (J53)l we arrive at 



^- + 777/' + 2r/ 2 = 2 (77^ + 2j/Vb) - vW - Zdr]' 
If we multiply the previous equation by 77, then add and subtract r/'r]" /2 we arrive at: 



dx v ' 1 2 2 rfx w ' ' ' ' " 2 

Integrating with respect to x and moving all the terms to the right hand side we get: 

// /2 4 

+ V 2 V' + \ - 2V ri 2 + drj 2 + c = (55) 

being c a real constant. Therefore, given Vq(x) the new potential V^x) and the function 
7(x) are obtained through (jlSj) and (|52~|) once we know explicitly a solution r?(x) of To 
get it, we make the Ansatz [65,66,77] 

77' = -r/ 2 + 2/??7 + 2£ (56) 

where (5{x) and £(x) are functions to be determined. Therefore: 

v " = -2r}T]' + 2/3't? + 2/V + 2f (57) 

= -77V + /?'?7 2 + /377T7' + 77^' (58) 

= -^ + ^ 3 +(e-/? 2 )r / 2 -2^-e 2 (59) 
The substitution of (|HH)l and (f5T)J) into (JH5|) leads to: 

p'r] 2 + (3t]7]' + 7?£' + /3t7 3 + (£ - /3 2 )r? 2 
-2/3^ - £ 2 - 2V ^ 2 + dr? 2 + c = (60) 
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and using again to eliminate rf in the previous equation we arrive at: 



{(3' + (3 2 - 2V + £ + d) t] 2 + - + c = (61) 

Since (16 lj) must be valid for arbitrary rj, the coefficients of each power of rj must vanish, 
which leads to £ 2 = c. Therefore: 

P'(x)+l3 2 (x) = 2[V (x)-e} e = ^ (62) 

Alternatively, we can work with the Schrodinger equation related to (jH^j) through the change 
(3 = [97]: 

_L(o)" + VqU (o) = m (o) (63) 

Depending on whether c is zero or not, £ vanishes or takes two different values £ = ±^/c. If 
c = 0, we need to solve one equation of form (}62|) and then the resulting ()56)1 for rj(x). If c 7^ 
there will be two different equations of type (|62jl. with factorization energies e± = (d + y/c)/2 
and e 2 = (d — y / c)/2. Once we solve them, it is possible to construct algebraically a common 
solution rj(x) of the corresponding pair of equations (|56|). There is an obvious difference 
between the real case with c > and the complex case with c < 0; thus, it follows a natural 
scheme of classification for the solutions r)(x) based on the sign of c. 



3.2 Classification of the second-order SUSY transformations 
3.2.1 The real case with c > 

Here we have that ei,e 2 G K, e\ 7^ e 2 , and the corresponding Riccati solutions of (J52"j) are 
denoted by Pi(x), $2 (a;) respectively. The associated pair of equations (J5fi)l become 

t/( x ) = -rf(x) + 2/? 1 (x)r/(x) + 2(ei - e 2 ) 
t/(x) = -?f(a;) + 2p 2 (x)r]{x) + 2(e 2 - ei) 

By subtracting them we get algebraically r/(x) in terms of e 1; e 2 and Pi(x), f3 2 (x): 

**> ^ "Jot ( 64 > 

If the corresponding Schrodinger solutions are used we have 

H) 4 0) ) w(»f , ,4 0) ) 

where W(f, g) = fg' — gf denotes the Wronskian of / and g. It follows that the second-order 
SUSY partner potentials V 2 (x) have no added singularities if W{uf' , u^) has no zeros. 
The spectrum of H 2 , Sp(H 2 ), will differ from Sp(-ffo) depending on the normalizability 

(2) (2) 

of the two mathematical eigenfunctions ipei , V'q of H 2 associated to e\ and e 2 which belong 
to the kernel of B: 

B4f = H*l>W = e j 1>® J = 1,2 
13 



For the solution associated to t\ we explicitly have 

0' + V0 + (7 + v'Wt? = (66) 
4f = 2(V 2 - e x )^ (67) 

By substituting (ffiTj) in one finds: 

mPif + + V + 2V 2 - 2e 1 )if>V> = (68) 
and using the expressions for V 2 and 7 given in (J48j) and (|52j) with d = €\ + t 2 we get: 

4? r,'/2 - V 2 /2 + e x - e 2 



(69) 



But from our Ansatz 
so that 

Therefore 

A similar procedure leads to 



rf = -rj + 2/3xri + 2{ex - e 2 ) (70) 

;(2)' / / (0)' 

(0) 

(72) 

,( 2 ) 



W(r4 0, ,t4 0) ) 

Concerning the possibilities for manipulating spectra offered by the second-order super- 
symmetric quantum mechanics, we have found a heuristic criterion providing some interesting 
information. Remember that 

B 2 Bl = (H -e 1 )(H -e 2 ) 

which has to be positive definite on 7i. In particular, this has to be valid on the basis of 
energy eigenstates \tp®) of H , and thus we have 

MB 2 Bt\<4>° n ) = (E n - ex)(E n -£ 2 )>0Vn 

This opens unexpected possibilities for the positions of the new levels e 1; e 2 . A non exhaustive 
list of several interesting situations useful for the spectral design is presented next. 

a) Our heuristic criterion indicates that if e 2 < e x < E it is possible to find «j° and 
such that W(u\^\u 2 ^) is nodeless and ipif are normalizable. Indeed, with 
this ordering of e\ and e 2 the right choice is to take nodeless and having one 
zero at xq, i.e., uf\xo) = 0. Since W'(u^\u2 ) = 2(ei — e 2 )u^ 'u 2 °\ it turns out 
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Figure 1: Second-order SUSY partner potential V 2 (x) (black curve) of the harmonic oscillator 
(gray curve) generated by using two Schrodinger solutions u 2 , of the form (|H7|) with 
e 2 = —1.2, v 2 = 1.1 and e x = —1, ^ = 0.9 respectively. Here Sp(if 2 ) = { — 1.2, —1, |, |, . . . } 



that W{uf\u { 2 0) ) has one critical point at Xq. Moreover, we have that at this point 
W"(uf \ u 2 ^) /W(uf\ uf^) = 2(ei — e 2 ) > which implies that the Wronskian acquires 
either a minimum positive or a maximum negative value at xq. Hence, W(%i[\ u 2 °^) 
is nodeless. The spectrum of the new Hamiltonian is Sp(if 2 ) = {e 2 ,€i, E , E 1 , . . .} 
(see a potential V 2 (x) with this kind of spectrum in figure 1). Notice that this case 
coincides with the one typically discussed when the second-order transformation is 
achieved through the iteration of two first-order transformations. 

b) If we choose Ei < e 2 < ei < E i+ i, our heuristic criterion suggests that perhaps we 
can find u± and u 2 such that W(uf\ u 2 ^) is nodeless and , ipi? are normalizable 
[96,98-100]. We would get then that Sp(H 2 ) = {E , . . . ,Ei, e 2 , ei, E i+i , . . . } (see figure 
2). This possibility becomes true if we select u 2 and having i + 2 and i + 1 nodes 
respectively. Indeed, due to the oscillatory theorem, which can be translated as the 
fact that between two zeros of there is at least one zero of uf\ it turns out that 
these 2i + 3 nodes are alternating. These zeros, ordered as x < X\ < ■ ■ ■ < x 2i+2 , are 
as well the critical points of W{uf\ u^) = W{x). Since W(xj)/W(xj+i) > 0, then 
W(x) conserves its sign (i.e., it is nodeless) in the interval (xj,Xj + i). Therefore, it 
does not have zeros in the domain (xo, £ 2 j +2 ). Finally, since xq is a node of u\ , then 
W"(xq)/W(xq) = 2(ei — e 2 ) > 0. Hence, W(x) acquires either a maximum negative or 
a minimum positive value at xq. In both cases W(x) never crosses the x-axis in the 
interval (— oo,Xq). A similar treatment leads us to conclude that the Wronskian does 
not vanish in (x 2 j +2 , oo), and therefore it is nodeless in the full real line. 

c) If e 2 = Ei, e x = E i+1 , and uf = 4°\ uf ] = 4% then W(uf\uf) is nodeless but 
tpei\ are non- normalizable. In order to prove this, let us notice that u 2 and 

have now i and i + 1 nodes respectively. Due to the null asymptotic behavior at 
x — > ±oo of both uf\ u 2 °\ it turns out that these 2i + 1 zeros, ordered as x < X\ < 
■ ■ ■ < %2i, are once again alternating, with xq and x 2 i being now zeros of . Using a 
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Figure 2: Second-order SUSY partner potential V^x) (black curve) of the oscillator (gray 
curve) generated through two solutions u 2 °\ of the form (JHTj) with e 2 = 3, z/ 2 = 0.9 and 
e x = 3.2, ui = 1.1 respectively. We have that Sp(H 2 ) = {| , §, §, 3, 3.2, |, . . . }. 

reasoning similar to that used in the previous case, W(x) becomes nodeless in (x , X2i). 
On the other hand, W(x) is monotonic in the interval (— oo,xo) and W"{xq)/W{xq) = 
2(e 2 — ei) < 0, which implies that W^x) has either a maximum positive or a minimum 
negative value at x = Xq. Since lim x _ t _ 00 W(x) = 0, it turns out that the only node 
which W(x) has in the interval (— oo, x ) is an asymptotic zero as x — > — oo. A 
similar procedure shows that in (x 2 j,oo) PU(x) has an asymptotic zero as x — > oo. 
In conclusion, the Wronskian is nodeless in the full real line, except by asymptotic 
null behaviors at x — > ±oo. This implies that the second-order SUSY transformation 
is non-singular inside the initial domain of definition. As the intertwining operator 
respect the boundary conditions of the eigenfunctions of H (except that now the 
eigenf unctions of H 2 associated to E{ and Ei + % are not square-integrable anymore), it 
turns out that Sp(iJ 2 ) = {-^o, • • • , -Ei-i, E i+2 , . . . }, i.e., somehow we have 'deleted' the 
two levels E i+ i in order to generate V 2 (x). 

According to the standard SUSY treatment, for which the new levels are always below 
the ground state energy of the initial Hamiltonian, the previous cases (b) and (c) are some- 
how unexpected, supplying us with more freedom for manipulating spectra. In principle 
these non-typical cases can be achieved through first-order SUSY transformations, but the 
corresponding interpretation is strange: in the first step we generate a singular potential 
V±, with the corresponding singularities induced by the zeros of the transformation function 
which is employed. The second transformation removes then the singularities introduced in 
the first step to arrive at a final non-singular potential V 2 . Next we will explore another 
interesting cases having unexpected positions for the new levels. 

3.2.2 The confluent case with c = 

In this case £ = 0, therefore e = e\ = e 2 G M.. Once we have found a Riccati solution /3(x) to 
for the given e, we must solve the Bernoulli equation resulting from (jSUjl [69, 101] 

7]' = —rj 2 + 2(3rj 
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In order to solve it, let us take r] — 1/y. Then 



V +2f3y = l 

Hence 

wq being a real constant. Thus, the general //-solution is given by 

g 2 / /3(x)dx 



r](x) 



w + J ePSPW^dx 
In terms of the Schrodinger solution u^(x) oc exp[J (3(x) dx] we have 

ju^(x)] 2 

V{ ] w Q + [: i y%y)fdy w(x) 
where xq is a fixed point in the domain of Vq(x) and, up to an unimportant constant factor: 

w(x)=w + [ X [u^(y)] 2 dy (73) 



In order that ^(x) will not have singularities, w(x) must be nodeless. Since w(x) is a non- 
decreasing monotonic function, a simple choice [101] (see also [95]) is to use Schrodinger 
solutions such that 

POO 

lim u {0) (x) = and 1+ = / [u (0) (y)f dy < oo (74) 

or 

/x 
[u {0) (y)] 2 dy < oo (75) 
-oo 

In both cases it is possible to find a ifo-domain where w(x) is nodeless [101]. For instance, 
if ({71} is valid and it^x) is a non-physical eigenfunction of if associated to e, it turns out 
that \im x ^_ oc w(x) = — oo and lirn a ._ j . 0O = w + J + . Thus, the domain for which 
is nodeless is < A similar procedure implies that for non-physical transformation 

functions satisfying (|73j) the nodeless w -domain is w > I_. 

(2) 

Similarly as in the case with c > 0, it can be found now a function ip\ in the kernel of 
B which is simultaneously an eigenfunction of H2 with eigenvalue e. Indeed, let us notice 
that equations (j6(iH71j) remain valid in this case, we just have to substitute e\ and €2 by e, 
u± by and (3i by (3. Thus: 



l(2) , . T}(x) U (0) (x) 



m(°)(x) w(x) 

(2) 

The spectrum of H2 depends on the normalizability of ijjl ■ In particular, for e > E$ it 
is possible to find solutions satisfying (f?4|) or (fT5|) and such that ipf^ 1 is normalizable. 
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Figure 3: Confluent second-order SUSY partner potential Vzfe) (black curve) of the oscillator 
(gray curve) generated by using a solution u^°' of the form (J87j) with e = 4, v = 1 and taking 
w = 5, x = in (73J). Now Sp(# 2 ) ={§,§, f , |, 4, |, . . . }. 




Figure 4: Confluent isospectral second-order SUSY partner potential V 2 (x) (black curve) 
of the oscillator (gray curve) generated by using the first excited state = ^ of the 
oscillator, with e = Ex = §, Wq — 1 and x Q = 0. 

This means that the confluent second-order SUSY transformations allow the embedding 
of single energy levels above the ground state energy of H (see figure 3). In addition, 
the physical solutions associated to the excited state levels of H can be used as well as 
transformation functions, no matter the number of zeros they have in the domain of Vq{x) 
(see figure 4). These two features cannot be achieved by iterating the first-order SUSY 
transformations without paying the price of introducing singularities at the intermediate 
potentials. As expected, these atypical spectral possibilities are consistent with our heuristic 
criterion previously formulated. 

3.2.3 The complex case with c < 

If c < we have that t\ G C and e 2 = e\. Notice that our heuristic criterion allows this 
possibility without violating the positive nature of the operator B 2 B\. We are going to 
analize just the case for which V^x) is real valued, implying that {5%{x) = /3i(x). Following 
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Figure 5: Real isospectral second-order SUSY partner potential V^x) (black curve) of the 
oscillator (gray curve) generated by using a solution of the form (J87j) with V\ — — 1 and 
ei = 5.4 + z/20 G C. 



analogous steps as for the real case, one arrives to the solution r\{x) of in terms of the 
(complex) solution $\{x) of the Riccati equation (|62)) associated to ei [102] (for the case 
when V" 2 (x) is allowed to be complex see [103,104]): 

v(x) = _ 2Im ^i) ( 76 ) 
Using the corresponding complex Schrodinger solution u\ (x) we can write 

m^) = 7~T- W(x) = : r— 

/l ' w(x) 1 ; 2( £l - h) 

In order to avoid the creation of singularities in V^x), w(x) must be nodeless. Since 
w'(x) = |u^(s)| 2 , it turns out that w(a;) is a non- decreasing monotonic function. Thus, to 
ensure that w(x) ^ V x G R it is sufficient that [102] 

lim (x) = or lim (x) = (77) 

x— >oo x— >— oo 

For transformation functions obeying the condition ()77|) it turns out that V^x) is a real 
potential isospectral to Vq(x) (see figure 5). Let us notice that, as in the previous cases, 
this transformation could be achieved by iterating the first-order SUSY, but the interme- 
diate potentials would be complex. This is, perhaps, the reason why this case was almost 
unexplored in the past. From our viewpoint, however, it represents a very promising line of 
research for the next years (see e.g. [105]). 
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4 An example: the harmonic oscillator 



We are going to apply the previous techniques to the harmonic oscillator. Before doing that, 
however, it is convenient to discuss some generalities of a subset of nonlinear deformations 
of the Heisenberg algebra. We will realize later on the importance of such structures for the 
SUSY partners of the oscillator. 

4.1 Polynomial deformations of the Heisenberg algebra 

First of all let us remember the standard oscillator algebra 

[H, a*] = a) [H, a] = —a 
[a, a)} = 1 

for which the number operator is a linear function of H: 

N = a)a = H-\ 

On the other hand, the polynomial Heisenberg algebras of m-th order are deformations 
of the previous structure, where there are two standard commutation relationships 

[H, L ± ] = ±L± (78) 

and an atypical one characterizing the deformation: 

[L~, L + ] = N(H + 1) - N(H) = P m (H) (79) 

where a generalization of the standard number operator is given by N(H) = L + L~ . The 
corresponding systems are described by Schrodinger Hamiltonians 

ff=-l£ + V(,) (80) 

being L differential ladder operators of order (m + l)-th, N(H) a polynomial of order 
(m + l)-th in H factorized as 

m+1 

N(H) = J] (H - Si) (81) 

i=l 

and P m (H) a m-th order polynomial in H. The algebra (j78H81j) generated by {H, L~,L + } 
provides information on the spectrum Sp(if) of H [67,79, 106, 107]. Indeed, let us consider 
the solution space of the (m + l)-th order differential equation (the kernel Kl- of L~): 

m+1 

L~ ijj = L + L~ i) = Y[ {H - Si) tp = (82) 

i=l 

Notice that Kl- is invariant under H: 

L'Hijj = (H + l)L~ip = Hip e K L - V ip e K L - 
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(b) 



Figure 6: Possible spectra for the Hamiltonians satisfying (|78 |1 — (|81 |) with s physical extremal 
states. In case (a) s infinite ladders are obtained by acting with L + . In case (b) there are 
s — 1 infinite ladders, and a finite one (the j-th one) is built out of ij}£. taking into account 

(S3). 



Thus, it is natural to select as the basis of Kl- the functions which are simultaneously 
eigenstates of H with eigenvalues Si 

H^ Si = S^ £i (83) 

These represent the extremal states for the m+1 mathematical ladders of spacing AE = 1 
starting from Si. If s of these states are physically meaningful, {ipet, i — lj • • • > s}, then by 
acting iteratively with L + onto them s physical energy ladders can be constructed (see figure 
6-a). 

It could happen [67] that for the ladder starting from Sj there is a I G N such that: 

(L + ) ,_1 Vfe,^0 (L + )'Vte, = (84) 

Then 

L-(L+)%=L-L + (L + ) l -^ £j 

'm+1 \ 
'm+1 



\l (£, + i-Si)\ (L+y- 1 ^ = o 



i=i 



This implies that another root of (|81)l must have the form Sk = Sj + I, k G {s + l,...,m + 
1}, j G {1, . . . , s}. Hence, Sp(if) will contain s — 1 infinite ladders and a finite one of length 
/, starting from Sj and ending at S 3 + I — 1 (see figure 6-b). 

We conclude that the spectrum of systems described by polynomial Heisenberg algebras 
of order m can have at most m + 1 infinite ladders. Notice that pairs of ladder operators for 
the harmonic oscillator satisfying (|78H81)1 with m > can be constructed simply by taking 
L~ = aP(H), L + = P(H)a Jf , where a,a^ are the annihilation and creation operators of the 
oscillator and P{H) is a real polynomial in H [108]. These deformations are reducible since 
for the same system we already have ladder operators a, obeying a much simpler algebra. 
Here we will be mainly interested in the search of irreducible deformed algebras. 
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4.2 SUSY partners of the oscillator 

Let us consider the harmonic oscillator potential 

rf 2 

Vo(x) = — 

The corresponding Hamiltonian has an equidistant spectrum with eigenfunctions and eigen- 
values given by: 



^\x) = J^—=H n (x)e-^ E n = n + ^, n = 0,l,... (85) 



where H n (x) are the Hermite polynomials. In order to implement the several SUSY trans- 
formations discussed previously, we will solve the Schrodinger equation (jfiHJ) for an arbitrary 
factorization energy e, namely: 

_I M (0)" + ^L(o) = m (o) 
2 2 

2 

Suppose that = e~^h(x). Therefore: 

h" - 2xti + (2e -l)h = 
By changing variables y = x 2 we arrive at: 

d 2 h . dh , „ 

V + (6_!,) ^- aft = (86) 

which is the confluent hypergeometric equation with a = (1 — 2e)/4 and 6 = 1/2. The general 
solution of (|86|) leads to the general solution we were looking for: 



u^(x) = e a 



, 1 „\ r(a + ±) / 1 3 , 



(87) 



where e G C and iFi (a, b; y) is the confluent hypergeometric (Kummer) function: 

lflMi ,)^£l>±^ (88) 
r («) ^ r(6 + n) n\ 

Notice that, for e < 1/2 and \v\ < 1 the solution given in ()87|) is nodeless. 

The explicit expressions for the SUSY partners of the oscillator can be calculated using 
(|87j) and, in general, they are too involved to be shown here. However, we will present the 
simplest derivation of a family of potentials isospectral to the oscillator [6] , which originally 
was derived by Abraham and Moses through the Gelfand-Levitan formalism [95]. 
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Figure 7: The Abraham-Moses family of potentials (|59*j) as a function of x and continuous 
values of the //-parameter in the domain [—0.95, 0]. 

4.2.1 The Abraham-Moses potentials 

Let us make a first-order SUSY transformation employing the general Schrodinger solution 
(J57j) associated to e = —1/2, namely: 

u (°)(ar) = e^[l + i/Erf (z)] 

where 

Erffar) = 4= / e _t3 dt 
V ^ 

is the well known error function. The corresponding superpotential is given by 

/ x 2z/e 2 

ai(g ' c)=a;+ 0F[l + ^Erf(x)] 

which is precisely the one obtained by Mielnik through his generalized factorization [6]. The 
first-order SUSY partner potentials are given by: 

™ = t - (v^Sifwi) - 1 (89) 

Up to a displacement in the energy origin, this is the Abraham-Moses family of potentials 
isospectral to the oscillator [95]. A plot of these potentials in terms of the variable x for 
continuous values of the parameter v in [—0.95, 0] is given in figure 7. Notice that for v = 
we recover the original oscillator potential. 



23 



Let us perform now a confluent second-order SUSY transformation employing the ground 
state eigenfunction, namely: 

« ,0) = < - J 0») 

For the confluent algorithm the key u>-function becomes (see equation ()73jl): 



w\x) 



wq + |Erf(x) 



where we have made Xq = 0. Taking now 1/wq — 2v, the new potentials read: 



2 / l\ I 1 I n — — 

rr , . x / «/ \ x / 2z/e 2 

Uo (X) = — ~ 



wj 2 I v^F[l + i/Erf (x)} 

These are precisely the already mentioned Abraham- Moses potentials (compare with ([89))) . 
Indeed, the procedure of derivation based in the confluent algorithm employing (j9~Uj) essen- 
tially coincides with the Abraham- Moses treatment [95]. Although in this case the first-order 
SUSY and the confluent second-order SUSY algorithms led to the same family of potentials, 
this was a fortunate coincidence due to the huge symmetry of the oscillator. In the gen- 
eral case, the potentials derived through the first-order SUSY will not coincide with those 
generated through second-order SUSY transformations. 



4.3 Non-linear algebra of Hk 

Suppose now that we have applied a k-th order SUSY transformation to V = x 2 /2, producing 
then the potential by creating k new energy levels. The spectrum of the end Hamiltonian 
Hk, intertwined with the harmonic oscillator through B' k , will be {ei,E n = n + 1/2, i = 
k, . . . , 1, n — 0, 1, . . . }, i.e., it contains a part isospectral to the oscillator plus k additional 
levels 6j, i = k,...,l placed by simplicity below E = 1/2. Hence, it is possible that 
polynomial Heisenberg algebras rule the H-SUSY partners of the oscillator. 

To analyze the algebraic structure characteristic of the Hamiltonians H^, let us look for 
ladder operators which connect the eigenstates associated to the levels E n . There is a natural 
construction for a pair of these operators [6,67,84,85], which is guessed from equation (|3l)|). 
its adjoint and the standard intertwining relationship involving the oscillator Hamiltonian 
H and its creation and annihilation operators a' , a: 

(H - l)a f = a j H {H + l)a = aH (91) 

The construction is composed of three stages (see figure 8): i) first we 'move' the eigenvectors 
|-0n) of Hk, represented previously by the wave functions tpn \ to the eigenvectors of 
the oscillator Hamiltonian H through the intertwining operator Bk- ii) Then, we move up 
(iV'n+i)) or down (\ipn-i)) on ^ ne ladder of H by using a) or a respectively, which will cause 
the effective 'motion' up or down on the ladder of Hk- iii) Finally, we get back to the ladder 
of Hk by acting B\ on l^n+i) or Wn-\l- Thus, the 'natural' ladder operators for Hk are: 

Lk = B\aB k L\ = B\a ] Bk k = 0, 1, . . . (92) 
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Figure 8: Diagram representing the action of the fc-th order intertwining operators B k , B\ 



and of the ladder operators a, at, L 



L\ for the Hamiltonians Hq and Hk 



For completeness, we have extended the intertwining relationship (J3HJ) to the case k = 
by assuming that Bq = B = I, I is the identity operator. The action of L k and L\ is 
drawn just onto the spectral points E n = n + 1/2, n = 0, 1, . . . because the eigenstates 
), i — 1, . . . , k} are annihilated by both L k and since they are annihilated by B k . 
Notice that 

H k L k = H k B\aB k = B k H aB k = B k a(H — l)B k 
= B\aB k {H k -l)=L k {H k -l) 

i.e., L k and L\ are differential ladder operators of order [2k + l)-th satisfying: 

[H k ,L k } = -L k [H k ,L\} = L\ (93) 

The operator N(H k ) = L k L k , which generalizes the standard number operator for the har- 
monic oscillator, is a polynomial in H k of (2k + l)-th order [79,81, 109]): 

k 

N(H k ) = L{L k = BialB k BiaB k = B{a^l[(H -6 l )aB k 

i=i 

= Bl(H -'£}l[(H -e i -l)B k 

^ ' i=i 

^ ' i=i 
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Thus, the operators L k , L k and H k generate a polynomial Heisenberg algebra of order 2k: 

[L k ,L% = N{H k + l)-N{H k ) (95) 

For consistency, when k = we should get the standard Heisenberg algebra because L = a 
and L — a) . This linear case is indeed recovered from our formulae: 

[H , L ] = -L [H , 4] = L\ [L , Lj] = / (96) 

The corresponding number operator becomes the standard linear expression in terms of the 
oscillator Hamiltonian: 

N(H Q ) = Ho-\ = N (97) 

On the other hand, when k = 1 and e\ is arbitrary we recover the quadratic Heisenberg 
algebra [82] (see also [79,81]): 

[L 1 ,L\] = (H 1 -e 1 )(3H 1 -e 1 ) (98) 

The number operator becomes now cubic in Hi: 

N(Hi) = (Hi - |) (Hi - ei) (Hi - e a - 1) (99) 

If k = 2 we will get a polynomial Heisenberg algebra of fourth order: 

[L 2 , 4] = (H 2 - ti){H 2 - e 2 ) [5if 2 2 - 3(ei + t 2 )H 2 + ei e 2 + l] (100) 

where N(H 2 ) is a 5-th order polynomial of H 2 : 

N(H 2 ) = (H 2 - |) (H 2 - ei) (H 2 - e 2 ) (if 2 - Cl - 1) (if 2 - e 2 - 1) (101) 

For general k, it arises a polynomial Heisenberg algebra of order 2k whose properties are 
characterized by the (2k + l)-th order polynomial N(H k ) of ()94j) . 

It should be clear now why the roots of are {1/2, e*, e^H-l, z = 1, . . . , k}: our H-SUSY 
partner Hamiltonians H k have precisely k + 1 physical extremal states associated to the + 1 
roots (eigenvalues) {1/2, e i: i = 1, . . . , k}. As the ladder starting from 1/2 is infinite, it does 
not imply any restriction on the remaining roots of the (2k + l)-th order polynomial N(H k ). 
However, as the ladders starting from £j are of length 1 (they end again at the initial energy 

the other k roots become €{ + 1, i — 1, . . . , k, as those appearing in (|9*1|). 

An interesting point to be addressed next, concerning the non-linear nature of the poly- 
nomial algebras (|92H95|) . is that they can be partially linearized [67,85,86]. 

4.4 Linearization of the non-linear algebra of 

As pointed out previously, the k eigenstates of H k , = 1, • • • , k, are isolated between 

themselves and from the \ipn) because L k \tp^) = = 4lV'e )- Hence, it seems natural 
to look for a linearization on the subspace spanned by {|V , n)? n = 0, 1, . . . }. The method 
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consists in modifying the ladder operators Lk and L\ of (j^2*j) to get an action similar to the 
resulting one when the Heisenberg algebra generators are applied to the appropriate energy 
eigenstates [67,85,86]. As for the sub-basis {\ip^),n = 0, 1, . . . } the commutator [L fc ,Lj,] is 
already diagonal (see equation ([95))). we propose a modification which will convert most of 
the diagonal elements of [Lk, L\\ to 1, namely: 

L L = Blf(N)aB k 4 = Bl<ff(N)B k (102) 

N = a) a being the standard number operator for the harmonic oscillator, f(x) a real function 
to be determined and the subscript L denoting linearization. We ask that [Ll,L^ l ] = I on 
the subspace Ti>\ spanned by {\i/;*),n = 1, 2, . . . }. Notice that we leave open the possibility 
that [L L , L\]\i/)q) = c\i[>q), c ^ 1 e R. It is straightforward to show that 

[L L ,L[M) = [g(n + l)-g(n)M) 

where 

g{n) = n [f(n - l)] 2 Utt (n - e, - \) (n - c, + \) (103) 
g(n + l)-g{n) = l n = l,2,... (104) 



The general solution of the finite difference equation (|104j) becomes: 

y(n) = n + w(n) (105) 
w(n) being periodic with period 1, w(n + 1) = w(n), n — 1, 2, . . . Hence: 



V n \Ai=l\ n - e i -2)( n - e i + 2) 

Since w(n) takes the same value for all n — 1, 2, . . . , it is important just w = w(l). Moreover, 
f(n — 1) should be real which implies that w > — 1. Finally, the ladder operators we were 
looking for read [67]: 



N + 1 + w 
(N + 1) nti(^-e,, + |)(iV - e t + 2 



Ll = bIJ — ■ — s — ^ ^ — aS fc (107) 



4 = ^R/ . iV + 1 + w 5fe (108) 

fe w (AT+i ) nt 1 (iv- ei + i)(iv- ei + |) 

Although their explicit forms are more involved than the ones for Lk and L k (compare (|92j) 
with (|107|108|) ). however and 4 act simpler on the energy eigenstates \ip^),n — 0, 1, . . . 
(except by the case with = which is separately discussed): 

LM) = (1 -Mv^+^l^ti) = Vn + w + l\^ k n+1 ) (109) 

[L L ,4]|^) = (l + w5 n o)l^) (HO) 
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This representation is independent of k, i.e., of the order of the interwining operator used to 
go from H to H k . The modified algebra here derived coincides with the 'distorted' Heisen- 
berg algebra originally introduced to linearize partially the second-order Heisenberg algebra 
characterizing the Abraham-Moses potentials, where w > —1 is the distortion parameter [85] 
(see also [67]). We have shown that this distorted algebra is common to all the H-SUSY 
partners of the oscillator. Moreover, a 'complete' linearization on 7i> (the subspace spanned 
by {\ip*),n = 0, 1, . . . }) can be achieved by taking w = to obtain precisely the Heisenberg 
algebra representation, namely: 

lM) = 4l^> = v^+T|^ + i) [l l AM) = M (in) 

If w = — 1 we get once again the standard Heisenberg algebra on 7Y>i but the state |^q) is 
annihilated by both and L L , i.e., it has been isolated by hand of the rest of eigenstates 
of Hk. This isolation property is natural for the other k eigenstates = 1, . . . , k. 

Notice that for k = the explicit expressions for L L , L\ do not coincide with a, a): 



iY+1 + W a L[ = a\ N+1 + W (112) 
N+l L V N+l K ' 

Now we get a distortion of the Heisenberg algebra representation which maps the operators 
a, into Ll, l) L , by changing the matrix elements of a and a) without affecting the diagonal 
elements of [a, a'] except the one associated to \iPq), which becomes w + 1. For w = we 
recover the original Heisenberg algebra because now Ll = a, L L — a). Moreover, when 
w = — 1 we get a reducible representation composed of the Heisenberg algebra on 7i>i and 
the null representation on the subspace generated by \i/jq) due to L l \iPq) = -kjrjV'o) = 0- 
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5 Coherent states for the SUSY partners of the oscil- 
lator 

The beautiful properties of the harmonic oscillator coherent states (CS) motivated the in- 
terest in looking for them in other physical situations [84-94,110-120]. There are various 
definitions, each one of them leading in general to different sets of CS. Concerning the inter- 
twining technique, CS which are eigenstates of a certain annihilation operator for the shape 
invariant potentials were constructed by Fukui and Aizawa [87]. The CS as eigenstates of the 
annihilation operator of ()92j) with k — 1 and e = —1/2, i.e. for the Abraham-Moses fam- 
ily of isospectral oscillator potentials were derived in 1994 [84], and the linearization process 
in the same case as well as the corresponding CS analysis was elaborated in [85,86]. Since 
then, a lot of works have arisen looking for interrelations between CS and quantum groups, 
pseudodifferential operators, non-linear algebras, etcetera [88-94,118,119]. In particular, 
the CS construction for the H-SUSY partners of the oscillator with arbitrary k has been 
successfully addressed [67]. For didactic purposes, let us discuss first the main properties of 
the standard CS [110-114]. 

5.1 Standard coherent states 

There are three equivalent definitions of the harmonic oscillator coherent states which can 
be used to define them in other physical situations. 

1. The coherent states \z) are eigenstates of the annihilation operator 

a\z) = z\z) zeC (113) 

2. The coherent states \z) are obtained by applying the displacement operator D(z) onto 
the oscillator ground state |0) (in this subsection we use the standard notation (Fock) 
for the eigenstates of the number operator, i.e., \n) = l^n))' 

\z) = D(z)\0) D{z) = e za "-~ za (114) 

3. The coherent states \z) are quantum states with a minimum-uncertainty relationship 

(Ax)(Ap) = \ (115) 

where for a system in the state \z) and an arbitrary operator O the corresponding 
uncertainty expresses: 

AO = ^/Wp - (z\0\z)) 2 \z) = ^{z\0 2 \z) - (z\0\z) 2 (116) 

Since for general systems we are specially interested in the first definition, we derive here 
the standard CS through ()113|) . First we expand \z) in the Fock basis 

oo 

\ z )=Y^c n \n) (117) 

ra=0 
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Then we use f!113|) to obtain: 

oo oo oo 

a\z) = y^c raA /n|n - 1) = y^zcjra) = zc n _i\n - 1) 

n=l n=0 n=l 

Hence, we get a recursion formula for c n : 



z Cn-l 



:ii8i 



n 



By iterating this equation we arrive at: 



Cn = ^ (H9) 



Finally, by asking that cq is a positive constant such that (z\z) = 1 we get: 

oo „ 
2 ^ 2, 

\z)=e- r -Y,^\n) (120) 



n=0 % 



where r = \z\. Let us remark that the standard coherent states can be used as an alternative 
basis (non-orthogonal) in the Hilbert space of states because they form a complete set: 

lf\z){z\#z = I (121) 

This can be straightforwardly verified since in the polar representation z = re td we have: 

- [ \z){z\d 2 z = - Y r e- r \ n+m+1 dr f 2 \ ie{n ~ m) d6 

*J *r t±n VM. Jo Jo 



n,m=0 



OO I \ i I poo 

^^p 2 y e - r \ 2n+1 dr = Y J W){n\ = I (122) 



n=0 " u n=0 



Thus, any state can be expanded in the basis of coherent states. In particular, any CS \z') 
admits a non-trivial decomposition: 

\z')= 1 -J\z)(z\z'}d 2 z (123) 

where the reproducing Kernel (z\z') is given by 

(z\z') = e -^ + ~ zz '- r -^ (124) 
Let us notice that coherent states evolve into coherent states, namely 



2 ytti 

itH \ z ) = e - r -yZ e -^H)t\ n ) = e~%{t)) z(t)=ze 



n=0 



30 



5.2 Coherent states for Hk 

Let us construct the CS for the SUSY partner Hamiltonians H k of the oscillator as eigenstates 
of the annihilation operators L k and L^. First, let us find them as eigenstates of L k (the 
non-linear case): 

L k \z) = z\z) zeC (125) 

We express \z) as a linear combination of the subset of eigenstates \ip%) °f H k associated to 
the part of the spectrum isospectral to the oscillator: 



^> = E c ^> ( 126 ) 



n=0 



After inserting ()126|) in ()125f) and using the fact that 



\ 



' ' 1\ / 1 



n 

8=1 



Win-ei--) n-e i+ - |V»*_i> (127) 



2y v 2 

we get a recurrence relationship for the coefficients c n 

v/ntir^ + Dr^ + DzVo 



n Ui=i( n - e i - \) \ n\Y[ k i=l T(n-e % + \)T{n-6i + 



:i28) 



We fix Co by the condition (z\z) = 1 and asking that c$ > 0. Hence, our non-linear CS 
become: 



' ' Z "" =0 y/n\aF 2k (-ei + ±,... ,-e fc +§ ,- £1 + |,...,- £fc + | ;r 2)n? =1 r(n- £i + |)r(n- £i + |) 1 J 

being r(a?) the gamma function, r = |z|, and pi 7 ^ a generalized hypergeometric function: 

rgi) ■ ■ ■ r(6 q ) ^ r(fll + w) . . . r( Qp + n) x" 

pF g (ax, . . . , a p , h, . . . , b q ; x) = —— — — > — — — r — — — - — - 130 

T(ai) . . . T(a p ) ^ T(bi + n) . . . T(b q + n) n\ 

Notice that z = is a (k + l)-th degenerate eigenvalue of L k because we get of ()129|) that 
\z = 0) = I^q) an d Lk\ip£.) = 0, i = 1, . . . , k. Thus, the resolution of the identity should be 
looked for as: 

k 

+ / (i3i) 



1=1 



where the measure d/j,(z) is to be determined. Suppose now that 

d/2(z) = F 2k + \, . . . , -e k + §, -ej + §, . . . , -e k + |; r 2 ) h(r 2 )rdrd6 (132) 

Inserting this equation in (j!31|) and using the fact that (iV'e-)) l^n); * = 1, n = 

0, 1, ... } is complete, we arrive at the following requirement for h(x): 

n h( w r(n + l)nt i r(n-e t + i)r(n-e t + f) 
x h(x)dx = t — — (133) 



^ntir(-6 i + |)r(-e i + 
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Hence, h(x) is the inverse Mellin transform of the right hand side of (j!33|) [121]. It turns 
out that h(x) is proportional to a Meijer G-function [67,121,122]: 

h ^ _ G?2fc o" 1 2fc+i( a; |0> ^finlllll! ~ ek ~ li ~ 61 + li • • • ! ~ efc + jj (134) 



^ntir(-e i + |)r(-6 i + 



Let us notice that for k = 1 and ei = —1/2 an explicit expression for ()134|) has been derived 
in [84]. 

The other properties of the standard coherent states have their corresponding analogue 
here. For instance, any CS of the form ()129j) can be expressed in terms of the others: 

\ z >) = / \ z )(z\z')dfi(z) (135) 



where the reproducing Kernel (z\z') can be straightforwardly evaluated: 

/~|~/\_ oF 2k (-e 1 + \,...-e k + \-t 1 +'± ,...,-e k + %;zz') (136) 

Vo F 2fe(-^l + i----^ + ^-ei + |v---^ + |;r2)oF 2fc (-e 1 + i,...-e fc + i -£! + §,... -e fc + |;r'2) 

This means that any two CS \z) and \z') of form ()129|) are non-orthogonal. From the 
resolution of the identity (|131|) it is clear that any state vector can be expanded in terms of 
our CS if we include the atypical orthogonal CS |?/^.), i = 1, . . . , k naturally inherent to this 
treatment. Notice also that our CS evolve in time as coherent states: 

U{t)\z) = e -%\z(t)) z{t) = ze- it 

Let us derive now the coherent states for the linearized annihilation operator Ll of (|107|) 



as: 



L L \z,w) = z\z,w) (137) 

where we have made explicit the CS dependence on the distortion parameter w. The expan- 
sion of \z, w) in the sub-basis |^) is equal to the one arising in the right hand side of (|126|) . 
and the use of (|137|) leads to the following recurrence relationship for the c n 's: 

ZCn-l y/T{w+l)z n 

+ w r(n + w + l) 

The coefficient Cq is determined by asking that (z,w\z,w) = 1 and Cq > 0. Finally 



\ Z}W) = r ( w+1 ) - Y - = =m ass) 

1,7 V ^(1, w + 1; r 2 ) ^ ^r(n + w + 1) ' Vn ' V ' 

The resolution of the identity becomes similar to ([131)1 : 

I = Y,\<^<\ + I \*,v)(*MdHL(z) (139) 
i=i J 
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where the measure dfii(z) is given by: 



duAz) — a(r,w) r dr d9 a(r,w) — - — — -^—4 — -e r r 2w (140) 
r w v ' ' y ' vrr(w + 1) v ' 

The reproducing Kernel is now 

/ 1 / v iFi(l,w+ l;g^) 

U,w z ,w) = . (141) 
ViFi^w+ljr^a.w+l;^) 

For w = (the linearized case in 7i>o) the same formulae as for the standard coherent states 
are recovered by noticing that 1 F 1 (l,l;r 2 ) = e r . Moreover, by taking carefully the limit 
w — > — 1 (the linearized case now in TC>i) it can be shown that the standard expression for 
the CS is also recovered, but the eigenstate l^o) associated to the eigenvalue E = 1/2 will 
be isolated of the other ones, i.e., the series ()138j) will start from |^) [67,85,86]. 

A comparison of the annihilation operators Lj~ and Ll and of both sets of coherent states 
derived in this section shows the following: the explicit expression for the non-linear operator 
Lk is simpler than the one for Ll- As can be seen from equations ()129|) and (j!38|) . however, 
the CS associated to L L are less involved than the ones associated to L k , which is due to the 
simplest algebra representation generated by L L and L L . 

Let us analyze, for both kinds of CS, the uncertainty products (Ax)(Ap) in cases for which 
we get closed analytic expressions. This can be directly achieved if we restrict ourselves 
to H-SUSY partner potentials in the harmonic oscillator limit. In particular, for k = 1 
and e± = —1/2 the uncertainty product for the non-linear CS ()129j) was calculated for the 
Abraham- Moses potentials in the oscillator limit [84], i.e., by taking v = in (j89|) . We shall 
present this case as well as the non-linear one with k = 2, e\ — — 1/2, e 2 = —3/2, and taking 
v\ = 0, i>2 — > 00 in order to recover the oscillator limit. For the linearized coherent states we 
shall study the case with k — 1, e\ — —1/2, v\ — 0, w = 1 (which was analyzed in [85,86]) 
as well as the one with k = 2, t\ = —1/2, e 2 = —3/2, v\ = 0, u 2 — > 00, w = 2. Notice that 
the corresponding expressions for the CS can be obtained from equations (|129|) and (|138j) 
by realizing that in the oscillator limit ~~ > lVn+2) anc ^ iVn) — y l^n+i)- Henceforth, the 
sets of coherent states we are dealing with read: 

\z) = , 1 E°°n f ) k = l (142) 



\ Z ) ~ yj F 4 (l,2 2 2,3;r2) E"=0 n!(n+1) t x /(^)i l^+2> & ~ 2 (143) 

u w= l) = i y« » _^ = uo ) fc = i (144) 

1 ' 1 ^/iFi(l,2;r2) ^"=0 ^/(n+1)! 1 r ™+ 1/ V ; 

l^ w = 2 )= v^^^Tfc 1 ^ fc = 2 (145) 

A direct calculation using ()142|) and x = (a + a))/y/2, p = i(a) — a)/y/2 leads to the 
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Figure 9: The uncertainty product (Ax)(Ap) (vertical axis) as function of z for the non-linear 
coherent states (jl42j) associated to the oscillator 



uncertainties Ax and Ap in the non-linear case for k = 1: 



Ax=\l--[Re{z)f Pl {r) 



Ap = 
Pi(r) 



- - [lm(z)] 2 Pl {r) 



0^2(2,2;^ 
0^2(1, 2;r 2 



0^2(2, 3; r 2 ) 
0*2(1, 2; r 2 ) 



(146) 
(147) 
(148) 



where Ke(z) and Im(z) represent the real and imaginary parts of z respectively. A plot of 
the corresponding uncertainty product is given in figure 9. 

On the other hand, for k = 2 and employing the non-linear CS in ()143j) it is found that: 



Ax = 
Ap = 
P2{r) 



- [Re(z)] 2 p 2 (r) 



- [lm(z)] 2 p 2 (r) 



F 4 (2,2,3,3;r 2 
F 4 (l,2,2,3;r 2 



oF 4 (2,3,3,4;r 2 ) 
6 Lo^4(l,2,2,3;r 2 ) 



(149) 
(150) 
(151) 



A plot of (Ax)(Ap) is given in figure 10. 

On the other hand, in the linearized case with k = 1 and making use of the CS in ()144j) 
we get: 



(Ax) 2 = (Ap) 2 = (Ax)(Ap) = - + 



2 iFi(l,2;r 2 ) 

A plot of the uncertainty product is given in figure 11. 

Finally, in the linearized case with k = 2 and employing ()145|) we find: 

1 



(152) 



(Ax) 2 = (Ap) 2 = (Ax)(Ap) = - + 



2 1^(1, 3; r 2 ) 



(153) 
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Figure 10: The uncertainty product (Ax)(Ap) (vertical axis) as function of z for the non- 
linear coherent states (j!43|) associated to the oscillator 




Figure 11: The uncertainty product (Ax)(Ap) (vertical axis) as function of z for the lin- 
earized coherent states (|144j) of the oscillator 




Figure 12: The uncertainty product (Ax)(Ap) (vertical axis) as function of z for the lin- 
earized coherent states ()145|) of the oscillator 
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A plot of (Ax)(Ap) is given in figure 12. 

As we can see, (Ax)(Ap) have more structure in the non-linear than in the linearized 
cases. Indeed, in the linear cases (Ax) (Ap) have radial symmetry, they differ of the standard 
result (Ax)(Ap) = 1/2 just in a neighborhood about the origin, and they quickly approach 
the standard value when \z\ — > oo (see figure 11 and 12). This does not happens for the 
non-linear CS for which the asymptotic values of (Ax)(Ap) are in general different from 1/2 
and depend on the direction in which we move out of z = (see figures 9 and 10). 

We conclude that the annihilation and creation operators for the oscillator H-SUSY 
partner Hamiltonians which are more similar to the standard oscillator ones become Ll and 
l) L . They mimic quite well the annihilation and creation operators of the oscillator and lead 
to the standard CS expression in the cases when the distortion parameter takes the two 
values w = and w = — 1. Moreover, L L and l) L become exactly equal to a and a) when 
k = w = 0, and the corresponding CS are precisely the standard ones for the oscillator. 
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6 Conclusions and outlook 



We have discussed the possibilities for designing quantum spectra offered by the supersym- 
metric quantum mechanics. We have seen that the standard iterative method, in which 
the first-order SUSY transformations are used to construct higher-order ones, could induce 
the wrong conclusion that the new Hamiltonians will have always the new levels below the 
ground state energy of the initial one. Here, it has been shown that the direct procedure 
frees us of that belief, allowing to create under certain restrictions some new levels above the 
initial ground state energy. It is important to underline also the complex SUSY transfor- 
mations which can be employed to generate non-hermitian Hamiltonians with either purely 
real spectra or with some finite number of complex levels [102-104]. We think that this 
interesting line of research is worth to be continued. 

When the SUSY techniques are applied to the harmonic oscillator, some interrelations 
immediately arise with several subjects of mathematical physics as non-linear deformations 
of Lie algebras and coherent states [67]. It can be established also connections between 
SUSY techniques and non-linear ordinary differential equations as Riccati [97], Painleve 
[75] and KdV [56]. In particular, when looking for the general systems ruled by second- 
order polynomial Heisenberg algebras, it turns out that the corresponding potentials become 
determined by a certain function which obeys the Painleve IV equation [107, 109]. It has been 
recently shown that subsets of H-SUSY partner Hamiltonians of the oscillator for which the 
k levels i — k, . . . , 1 are forced to be connected by appropriate annihilation and creation 
operators to form a finite ladder of equally spaced levels, supply us with explicit solutions 
of the Painleve IV equations [123] . Notice that a similar treatment relates subsets of SUSY 
partner potentials of the radial oscillator and solutions of the Painleve V equation [75]. 

We conclude by mentioning that links of this kind between SUSY techniques and non- 
linear aspects of mathematical-physics, as the ones previously pointed out, represent the 
future of a field which has proved very fruitful along the years (see e.g. [78]). 

Acknowledgements. The authors acknowledge the support of CONACYT (Mexico), 
project No. 40888-F. 
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